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The multi-flavor Schwinger model on at finite temperature T is mathe- 
matically equivalent to the model on 5^ at T = 0. The latter is reduced to a 
quantum mechanical system of — 1 degrees of freedom. Physics sensitively 
depends on the parameter m/T. Finite temperature behavior of the massive 
Schwinger model is quite different from that of the massless Schwinger model. 



1. Introduction 

The Schwinger model is QED in two dimensions, described by 
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Massless theory (m^ = 0) is exactly solvable. In the = 1 (one flavor) model the gauge 
boson acquires a mass without breaking the gauge invariance.^ It has the 9 vacuum 
and a non- vanishing chiral condensate (ipip) 7^ 0.^ 

The N > 2 model is distinctively different from the = 1 model. The spectrum 
contains — 1 massless bosons. The chiral condensate vanishes, (■?/''?/') = 0,^ as in 
two dimensions continuous symmetry, SU (N) chiral symmetry in this case, cannot be 

spontaneously broken. Nonvanishing (ijjip) in the A^=l theory is allowed, since the 
U{1) chiral symmetry is broken by an anomaly. 

With massive fermions the model is not exactly solvable. The effect of the fermion 
mass in the A^ = 1 theory is minor for m/e <^ 1, except that it necessiates the 6 
vacuum. For N > 2 the situation is quite different. Coleman showed"^ that in the 
N = 2 model with = m two resultant bosons acquire masses given by 
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Surprising is the fractional power dependence on m, e and | cos of the second boson 
mass /i2, resulting from the self-consistent re- alignment of the vacuum against fermion 
masses. The effect of fermion masses is always non-perturbative even if m/e <^ 1. 

""^To appear in the Proceedings of the 4th International Workshop on Thermal Field Theory and 
their Application, Dalian P.R. China, August 5 - 10, 1995. 



How does /i2 depend on m in the N flavor model? How does it change at finite 
temperature? Does the chiral condensate vanish at suflieciently high temperature? 
Can the fermion mass be treated as a small perturbation at high temperature? 

These are questions addressed in this article. We present a powerful method to 
evaluate various physical quantities at zero and finite temperature. We shall recognize 
the importance of a dimensionless parameter m/T. The behavior at m/T ^ 1 is 
quite different from that at m/T <^ 1. Coleman's result corresponds to m/T ^ 1 
as T ^ 0. At high temperature m/T <^ 1, one finds /i2 oc m, a result obtained in 
peturbation theory. This work is based on the result obtained in ref. 5. 

2. At finite temperature and on a circle 

In Matsubara's formalism the model at finite temperature T in equilibrium is equiv- 
alent to an Euclidean field theory, or a theory with an imaginary time r, satisfying 
boundary conditions 
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ija{r + -,X) = -iJa{T,x) , (t + - , x) = (t, x) . (3) 

If one, instead, places the model on a circle of circumference L (at zero temperature) 
with boundary conditions 

i:a{t,x + L) = -iJa{t,x) , A^{t,x + L)=A^{t,x) , (4) 

then one obtains a theory which is, after Wick's rotation, mathematically equivalent 
to the finite temperature field theory defined by (|^). Various physical quantities in the 
Schwinger model at T 7^ are related to corresponding ones in the model on by 
substitution of T by L^^. 

Our strategy is to solve the model on a circle with an arbitrary size L. There is 
powerful machinery which specifically works on S^.^~^ 

3. Reduction to a quantum mechanical system 

Fermion operators on 5*^ can be expressed in terms of bosonic operators: Take 
7^ = (cri,i(j2) and write ip'^ = {ip'^,'ip°L). In the interaction picture defined by free 
massless fermions 

^l{t,x) = ^C'^e±*^''i+2-Pi(*±-)/^> : e±^^<^±(*'^) : (5) 

The Klein factors are given by = e^''Sb=i and CI = e'^^^^^^^^P+'P-^ 
Here [g±,p±] = iS"'' and x) = E^=i(47™)-i/2 g-2™{t±x)/L + ^.c.} where 

[c±n5C±|m] = ^"'''Snm- The : : indicates normal ordering with respect to (c„,cjj). 
The antiperiodic boundary condition is ensured by a physical state condition (|). In 
physical states p^. takes integer eigenvalues. 
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After substituting the bosonization formula (^), the total Hamiltonian in the 
Schrodinger picture becomes 



Ho = + ^ E \iPl - P-f + iPl+P- + ^""'^ 
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Here p"^ — pt and p";^ + p"i correspond to the charge and chiral charge operators, 
respectively. Qw is the phase of the Wilson line around the circle, the only physical 
degreee of freedom of the gauge field on the circle, Ai = Qw{t)/eL. The coupling 
between p"^ and Qw is induced through the chiral anomaly/ 0a = 0+ + 0- and Ila is 
its canonical conjugate, ifmass is the fermion mass term. 

Notice that is an exact operator identity. In the absence of fermion masses 
Htot = Hq + H^. The zero modes {Qw, <1±) decouple from the oscilatory modes (pa, 
and the Hamiltonian is exactly solvable. The spectrum contains one massive field 
jY-i/2 Y^^^^ 0^ with a mass /i = (N/TcY^'^e, and — 1 massless fields. 

To examine effects of //mass, first note that -ffmass, and therefore -fftot, commutes 
with p";^ — p"l. Hence we can restrict ourselves to states with (p^ — p°L)\ phys) = 0. 
With this restriction a complete set of eigenfunctions and eigenvalues of Hq is 

r) N r) N 

^i"™^ = M^ + l) + f E-a-^{E-4 (8) 
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where a harmonic oscillator wave function Ug satisfies |(— 5^ + x^)us = {s + ^)us 
with X = {'Ke^L?' /N)~^/^Qw The ground states of Hq are infinitely degenerate for 
ni = ■ ■ ■ = riN due to the invariance under a large gauge transformation Qw OvF+27r 
and if) a — >■ e^'^^^^^ipa- 

Hrasss induces transitions among $("i'-'"^)'s. It also gives finite masses to the — 1 
previously massless fields. The structure of the vacuum sensitively depends on i^mass- 
The effect of ifmass turns out quite nonperturbative so long as nia ^ 0. 

It is more convenient to work in a coherent state basis given by 



{n,ra} 



Transitions in the s index may be ignored to a very good approximation. We seek the 
vacuum in the form 

$vac(^) = / dipi--- dipN-1 fi^a, 0) %{iPa; 0) . (10) 

Jo 



3 



Let Xa = Uaa4>a and fia be a mass eigenstate field and its mass. Then matrix 
elements of -ffmass in the coherent state basis are 

N-l N 
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B{z) is given by 
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The eigenvalue equation {Hq + ifmass) '^'vac(^^) = -E'$vac(^^) becomes 

{-A^ + \/;v(^)}/(¥^) = e/(^) 
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Here e = NEL/2tt{N — 1). Eq. (|T3|) is nothing but the Schrodinger equation with the 
kinetic and potential terms given by — and V)v(v^), respectively. 

The potential Vn{^p) depends, through Aa defined in (pUf), on fia and Uaa which are 
to be self-consistently determined from the ground state wave function f{ipa', 0) of the 
Schrodinger equation (0). /Xq's and UaaS are determined by 
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Ra = -j-Aa{ cos (fa) f = -4:7ima{i^ai'a) 9 ■ 



(14) 



We have denoted {g{'f))f = J[d<f] g{<f)\f{<f)\'^. We need to solve (|Tl|), (|T3D, and (|T|) 
self-consistently. 

We have shown that the N flavor massive Schwinger model is reduced to quantum 
mechanics of — 1 degrees of freedom in which the potential has to be fixed self- 
consistently with its ground state wave function. 
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Figure 1: Chiral condensate —{4'4')/fJ^ and as a function of temperature T at = 
in the = 1 and N = 3 models. In the N = 3 model m/ fj, — 10~^, 10~^, 10~^ and 

10"'^. The crossover in (-0-0)/// takes place when ~ 0.2. 



4. N=l (one flavor) 

One flavor case is special, (fi = 9 and there is only one massive boson with a mass 
Hi . This case was analysed in detail in refs. 7 and 9. 

The vacuum on is $vac(^) = ^o(^)- Converting all expressions to finite temper- 
ature case, we find the boson mass and chiral condensate to be 

= + SvrmT cos^ e"'^^/^ 
{^i;)o^-2Tcosee-^^/^B{^) (15) 
where = e/^/^T. At T = 

/ii = Y^//^ + m^e^'T cos^ 6 + me' cos 9 

- el 

{^^)^ = -—H^cos9 . (16) 
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Notice that the condensate is nonvanishing even for m = 0. As T is raised, it shows 
a crossover around T = /i, and approaches zero at high temperature. See fig. 1. The 
correction due to the fermion mass m <^ /i is minor. 

5. N>2 (multi-flavors) at low and high T 

The situation in the multi- flavor case is quite different. When rua = 0, Vn{(p) = 

and f{(p) =const so that ( cos ipa)f = 0. Consequently {ipa4'a)9 = 0. 

Suppose that fermion masses are degenerate: = m. There results one heavy 
boson and N — 1 light bosons with masses /xi and /i2, respectively. The potential in 
(|1^) becomes 



N 
a=l 



where /ii and /i2 are determined by 

2 2 2 

/i^ = + /i2 

f^2 = ——^ — — KoT"^ { cosLp)f = -4:Trm{^pai'a)e ■ (18) 

Recognize that two parameters kq and 6 fix the potential Vat (</)). 

If m 7^ 0, kq becomes very large at low temperature T —>■ 0. In this regime the 
potential term dominates over the kinetic energy term in the Schrodinger equation 
(|l^). The wave function f{ip) has a sharp peak around the minimum of the potential. 
The minimum is located at (p^ = d/N where = 0- 2-n[{0/2Ti) + |]. As varies from 
— TT to +7r, the minimum moves from = —n/N to ipa = +tt/N, and jumps back to 

ifa = -7t/N. 

In the T ^ limit, ( cos )/ = cos so that 

= -_f2e^cos^)'"+^ f-)^+^ for T < m^^/i— . (19) 
/i 47r^ N ^ n ^ 

Two important observations follow. First, the dependence of the condensate on m 
is non-analytic. It has fractional power dependence. The effect of fermion masses is 
nonperturbative in this limit. Secondly, as a function of 0, the condensate has a cusp at 
= ±7r, which originates from the discontinuous jump in the location of the minimum 
of the potential. 

In the opposite limit kq -C 1, which includes both m — (with T > kept fixed) 
and T 3> /i, the potential Vn{p) can be treated as a small peturbation in (^). One 
finds 

/ \ / (1 + cos6')ko for = 2 , . 

^''''^^^f = [ no for AT > 3. ^'^^ 
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Figure 2: Wave function \f{ip)\'^ in the = 3 model. 



There appears no 6'- dependence for iV > 3 to this order. The condensate for > 3 is 
found to be 



^'^^ 7r(iV- 1) /i 



g-2^T/7VM for T > 



For N = 2, the expressions for {%1)'iI))q in (|2l| ) must be muhiphed by a factor 2 cos^ 
Notice that the condensate is hnearly proportional to m in this regime. 

6. Crossover in {i'4')gT 

For general values of T and m, we need to solve the set of equations numerically. 
In the two flavor case the Schrodinger equation is 

f d'^ 1 

2ko cos ^9 cos(^i - m fiifii) = t /(^i) . (22) 



This is the equation for a quantum pendulum. The strength of the potential is 
2^0 cos which changes the sign at = ±7r. In particular, the potential vanishes 
at = ±7r and /(v^i) =constant. Consequently the chiral condensate vanishes at 
6'=±7r. This is a special feature of = 2. 
For = 3 the equation is 



COSV9i + COSV92 + COs(v9i + V92-^) jf{f)=^f{f) ■ (23) 

7 




Figure 3: The location of the minimum of the potential in the N = 3 model is displayed 
with various values of {q,r) in (23). For (g, r) = (10, 1.3) the minimum at 6'=±7r is 
located at the origin. 



The potential term never vanishes unless kq = 0, or equivalently m = or T ^ oo. 
The equation can be solved numerically for an arbitrary kq. The ground state wave 
function \ f{ip)\'^ has been displayed for various and 6 in fig. 2. 

With given values of m and T the chiral condensate is determined by solving (p!7D, 
and (|23|) simultaneously. We developed an iteration procedure which yields a 
consistent set of values of m//i, T/fi, fia/ fJ', and hq. It takes less than ten iterations 
even for moderate values of ~ 1 to achieve four digits accuracy. 

In the top figure of fig. 1 we have displayed the condensate {'ipip )e=o/ as a function 
of T/fi and m/fi in the A^=3 case too. With a given m/fi the condensate is almost 
constant at low T/fi, and sharply drops to a small value around T*. This crossover 
takes place when nQ{m^T^) ~ 0.2 for a wide range of the value of m//i. In the bottom 
figure we have shown a plot for in the same region of the m-T space. 

The asymptotic formulas (^) and (|2l|) are quite accurate for kq> 1 and < 0.1. 
The important parameter is Kq. 



7. The singularity at = ivr and fermion masses 



As shown in (|T^), the chiral condensate at T = shows a cusp singularity in its 
9 dependence when fermion masses are degenerate. When fermion masses are not 
degenerate, the potential Visfi'^) in (|13D is deformed from the symmetric one in (|1^. 
Accordingly the location of the minimum is shifted. The cusp singularity appears when 
the location of the minimum of the potential makes a discontinuous jump. 

With given {ma}, Aa in Vn{^p) is determined self-consistently. We have determined 
the location of the minimum of the potential in the three flavor case. The potential is 
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proportional to 

F{(pi, (p2;0) = -q cos (fi - r cos (p2 - cos(6' - (fi - (P2) (24) 

where q = A1/A3 and r = A2/A3. 

In the symmetric case q=r=l the location of the minimum moves from (— |7i", ~|7r) 
to (+^7r, +^vr) as 9 varies from — vr to +7i, and makes a jump. At ^=0 the minimum is 
located at the origin for arbitrary (g, r). We have plotted trajectories of the location 
of the minimum for several typical values of (g, r) in fig. 3. 

So long as asymmetry is small, there is a discontinuous jump ai 6 = ±7r. However, 
a sufficiently large asymmetry restores continuity. For instance, with {q < 0.5, r = 1) 
the minimum at 6' = ±7r is located at ^92) = (±7r, 0). The trajectory makes a 
closed loop on the (pi-(p2 torus. For (g ^ l,r = 1), the minimum at 6^ = tt is around 
(0, ^Tr) so that the discontinuity remains. However, if one adds a small aysmmetry in 
"r", the minimum is pushed back to the origin. For instance, for (g, r) = (10,1.3), 
the minimum starts to turn back at ~ O.Stt and reaches the origin at 6 = n. The 
implication to QCD physics is profound. As ^ ma > rriu, we are facing at a case 
g 3> r > 1 in which there is no singularity in 6 any more.^° 

We conclude that a sufficiently large asymmetry in the fermion masses removes the 

cusp singularity at6' = ±7rin {ipip) at T = 0. 
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